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Abstract 

The fluid model for the dark sector of the universe (darkon fluid), 
introduced previously in [T], is reformulated as a modified model in- 
volving only variables from physical phase space. The Lagrangian of 
the model does not possess a free particle limit and hence the parti- 
cles it describes, darkons, exist only as a self-gravitating fluid. This 
darkon fluid presents a dynamical realisation of the zero-mass Galilean 
algebra extended by anisotropic dilational symmetry with dynamical 
exponent z = |. The model possesses cosmologically relevant solu- 
tions which are identical to those of [1]. We derive also the equations 
for the cosmological perturbations at early times and determine their 
solutions. In addition, we discuss also some implications of adding 
higher spatial-derivative terms. 
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1 Introduction 



In the previous paper pQ we have presented a model of dark energy based 
on new particles which, in the following, we call "darkons" . These particles, 
which are nonrelativistic and massless, possess a modified relation between 
energy and momentum (or velocity) and so in their description we were forced 
to use an enlarged phase space (see [2]). The gravitational coupling of the 
darkons was introduced in a minimal way which satisfied the general form of 
the Einstein equivalence principle. This minimal coupling lead to a dynami- 
cally generated active gravitational mass density of either sign as a source of 
the gravitational field. Such a property opened the possibility of using this 
model to explain the observed accelerated expansion of the universe (see [3] 
for a very recent review). Our paper [1] then used a fluid mechanical gen- 
eralisation of this particle concept to construct a new parameter free model 
for the dark sector of the universe. And, when we performed a comparison 
of the predictions of this model with the cosmological observations, we were 
pleased to see no disagreement thus suggesting that we should look at it more 
seriously. 

The enlarged phase space, which was used in the construction of this 
model, seems somewhat unphysical. Hence we have decided to see whether 
we can reformulate the model using only the conventional phase space. This 
paper presents such a reformulation. 

The second important property of our model is the anisotropic scale sym- 
metry defined by the following transformation properties of space and time 
coordinates: 

f* = Ax, t* = X z t, (1) 

where z is a dynamical exponent. Any relativistic scale invariant theory 
clearly shows isotropic scaling (i.e. z — 1) as well as possesses the corre- 
sponding nonrelativistic limit |4j. Galilei covariant theories with nonvanish- 
ing mass can exhibit the anisotropic scaling with z = 2 [5], the so-called 
Schrodinger symmetry. Galilei covariant theories with z ^ 2 have necessarily 
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the vanishing mass [2] and so this is the case also for our model, which as 
we will demonstrate, corresponds to z — §. Such a z- value for a gravita- 
tional theory is exceptional. Recent renormalizable gravitational theories, 
which violate Poincare symmetry in the ultraviolet limit and which exhibit 
anisotropic scale symmetry there, usually correspond to z = 3 [6]. Cosmolo- 
gies with other z- values have also recently been discussed (see [7]). 

But there is a crucial difference between Horava's gravity [6J and our 
work in the way they get z ^ 1. In [6 J higher-order spatial derivative terms 
corresponding to z = 3 are introduced from the very beginning - in order to 
enforce renormalizability of the theory. In our case, however, the anisotropic 
scaling with z = | is an emerging symmetry resulting from the minimal 
coupling of massless nonrelativistic particles to gravity. 

For any model, to be a candidate for the description of the universe, one 
should determine, in addition to its cosmological solution, also the perturba- 
tions around this solution. This is what we do for our model in this paper. 
Basing our discussion on the symmetry properties of our model we derive 
perturbations around its cosmological solution at early times and present a 
preliminary discussion of their form. 

The paper is organised as follows. In Sect. 2 we review the previous 
formulation of our model and describe the modifications in its description so 
that the unphysical degrees of freedom used in its original formulation are 
absent. In Sect. 3 we show that the darkon fluid is an interacting system 
and, as such, possesses zero mass Galilean symmetry which will be extended 
to anisotropic scaling in Sect. 4. In Sect. 5 we derive the equations for cos- 
mological perturbations at early times and describe their solutions. Finally, 
in Sect. 6, we add higher spatial-derivative terms to the interaction and dis- 
cuss their implications. We conclude with some remarks and an outlook for 
further research in this area (Sec. 7). 
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2 The reformulation of the model 



2.1 Previous formulation 

The two-component fluid model in pQ is denned in terms of a Lagrangian 
given by 

L = L M + L D + L#, (2) 

where M, resp. D, stands for the baryonic resp. dark sector. 

In pQ we have started with another interpretation: M for the matter 
(baryonic and dark) and D for the dark energy sector. But as shown in 
|1] the D sector can be used for a unified description of the whole dark 
sector. Hence, in the present paper, we will follow this attitude from the 
very beginning. 

For the separate parts of L in fl2]) we have 

Lm = run™ J (V(£f - lyf) - , (3) 

where m is a mass parameter giving ([3]) the correct dimension, 

L D = n° J d 3 £ ( Pi (xf - yf) + qfijf + qf d t <P(x D ,t)) (4) 



and 



J d 3 x (V<p(x,t)y , (5) 



87TG 

where no denotes the constant particle distribution in £ space (£ is a contin- 
uous particle label) in the respective sector and G is Newton's gravitational 
constant. In these expressions all phase space variables are functions of £ and 
t, ie x M = x M (£,t) etc. Note that the variable q D ((,,t) acts as a Lagrange 
multiplier field □• For more details see pp. 

From (j4j) we note that the phase space in the D sector is now 12-dimensional. 

The equations of motion (EOM), corresponding to L, are then given by 



^or the role of Lagrange multiplier constraints in gravity theory see [5]. 
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• M sector 




,M 



yf = -d l( p(x M ,t) 



(6) 



• D sector 





n 



(7) 



yP = -d l( f>(x D ,t) 
pf = q k d k di(f)(x D ,t) 



• (f) sector 




(8) 



The last term in (jHJ) represents a dynamically generated active gravitational 
mass density of either sign leading either to an attractive or a repulsive 
gravitational force. For this reason the D— sector may serve as a model for 
the dark sector of the universe (for more details see pQ). 

All this discussion has been given in the Lagrange formulation. However, 
for the fluids it is more convenient to express all quantities in the Eulerian 
formulation. 

In such the dynamics of a fluid is usually described by the fluid 

particle density field n(x,t), the velocity field u(x,t) and the components of 
the force fields. These fields are given by 





and 



Ui(x,t) = ?/i(f(f,t),t), (10) 

where yi(£,t) is an independent degree of freedom and the functions £(x,t) 
are the inverse of i). In our case [1] we have, in addition to the n and 
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u fields, also the momentum Pi(x, t) and pseudo-momentum qi(x, t) fields in 
the D sector which are given by Pi(£(x, t), t), resp. qi(£(x,t),t). 

The corresponding EOM are then given by (for more details see [T]) 

d t n A (x,t) + d k (n A u A )(x,t) = 0, (11) 

where A e (M, D), ie the continuity equations for the particle number den- 
sities n M and n D . 

The Poisson equation for the gravitational field is given by 

A<t>(x,t) = 4nG (p M + di(n D qd) , (12) 

where the mass density p M is defined by p M := mn M . 

Note that the last term in ( II 2p represents the dynamically generated 
active gravitational mass density of the dark fluid. 

As shown in [TJ we have, in addition, the following equations: 

Dfuf = -d^ (13) 
from the second equation in ([6]) and from the third equation in (jTJ) 

Dfuf = -84, 

where we have defined 

Df = d t + u?d t . (14) 

Note that if uf 1 and uf obey the same initial conditions (|T3^) shows that 
uf = uf' 1 = Ui ie (fT5|) becomes one universal equation valid for all fluid 
components. 

D t m = -di<j). (15) 
Finally, the second and fourth equations in (J7|) give 

Aft = ~Pi, DtPi = qkdid k (f). (16) 

Looking at fl 1 5 |l 1 6 [) we note that, in contrast to standard fluid mechanics, the 
two vector fields p(x, t) and u(x, t) are not parallel to each other. 
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2.2 Modifications of the dark sector 

All this was introduced and discussed in [Tj. However, let us observe that 
the gravitational field (p(x,t) only appears in Lp and in the expression 
for the gravitational force gi = —di<p. This suggests to introduce gi as the 
field variable in the Lagrangian restricted to the dark sector and its self- 
interaction. Then, instead of fl8]) we would havaj 

g % (x,t) = -AnGrio J d 3 £5 3 (x- x{£, t) ) Qi {£, t), (17) 

which can be rewritten as 

gi(x,t) = -47rGn{x,t) qi (i(x,t),t), (18) 

and so 

*(£*) = - A(-) *V ( 19 ) 



AnG V n 

With these modifications our EOM become 



and 



where 



and 



D t Ui(x,t) = Qi{x,t) (20) 

= ^b) %fc(f ' t))2 ' (21) 

D t = d t + u k (x,t)d k (22) 



«(£,£) = £(f,% = ^)- (23) 

To complete the set of our equations we have to include also the continuity 
equation for the particle density 

3 t n + d k (nu k ) = 0, (24) 



2 As only the D sector is considered we omit the superscript D here and in the following. 
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which we generate by introducing a Lagrange multiplier field 9. Then the 
complete Lagrangian for the dark sector, in the Eulerian formulation, takes 
the form 

d 3 x (giiDtUi - \ gi ) + 9{d t n + ^(«»n) J . (25) 



4ttG J V 2 
Note that varying the Lagrangian with respect to the darkon density n 
gives us the equation for the lagrangian multiplier field 9: 

D t 9 = 0. (26) 

Furthermore, (12~T]) is obtained from the Lagrangian (125]) in two steps. By 
varying Ui we obtain the EOM 

D t Qi = gkdiU k - gid k u k - ndi9. (27) 

Then we apply D t to ( 127]) and obtain, by using the other EOM, eq.( l2T|) . Note 
also that in our solution (|T8|) of the Poisson equation for the gravitational 
field 

d i9i = -iirGdiinqi) (28) 

we have neglected a possible contribution to which is a curl of another 
vector field A(x, t). The neglect of such term (ie putting A = 0) is not a 
serious defect of our modified model as the choice (1T81) leads to the most 
general radial symmetric solution of the Poisson equation (1281) and so to the 
same cosmology as discussed in pp. 

Let us summarise our results: We have successfully reached our main aim 
of reducing the phase space to the physical one while keeping the cosmological 
solutions derived in [1] . However, we have had to pay a prize for this as now 

• the gravitational force is not automatically given by a gradient of a 
potential, 

• the EOM (1271) breaks the symmetry with respect to arbitrary time- 
dependent translations: 

Xi ^x[ + Oi(t), (29) 
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as, under this symmetry, the fields A should transform as A — > A' with 

u'i(x',t) = Ui(x,t) + di(t), 

9i(x',t) = gi(x,t) + 'di(t), (30) 

n'{x',t) = n(x,t) and 9'{x',t) = 9(x,t). 

Transformations with cij(£) being at most linearly dependent on t be- 
long to the group of Galilei transformations, whereas higher power 
ones (with d^O) are gauge transformations which represent Einstein's 
equivalence principle. 

It turns out that both problems can be solved by a further reduction of 
the dimension of the phase space. To do this let us assume that the velocity 
field u is curl-free, i.e. 

u = Vu. (31) 

Note that the restriction fl3T|) does not change the cosmological solutions. 
Next we insert (|3~T|) into the Lagrangian f )25|) and obtain 

L = y d 3 x ^(dtu + ^duu) 2 ) - X -g 2 % + eidtn + diindiu))) . 

(32) 

Varying gi we obtain again the EOM (1201) but with fl3Tl) the left hand side 
of (TSUI) becomes the gradient of a scalar and so g is derivable from a scalar 
potential </>: 

9 = -V0, (33) 
where is given by the Bernoulli equation 

-<t> = dtu+ \{d k u) 2 . (34) 

Varying the velocity potential u we obtain, instead of ( 1271) the EOM 

d t {digi) + di(d k g k diu) + <9;(n<9;0) = 0, (35) 
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which is nothing else than the divergence of (j27p if we take (13 ip and ( I33p into 
consideration. The EOM (|24l) . ( 1261) . ( 1341 and (1351) now constitute a system 
of four equations for the fields n, 9, (j) and u. 

It is easy to see that ( )35l) is now invariant under the set of transforma- 
tions ( l2T)f3"0~p . Note that the other equations ( KUp . flU]) and (12"6"P possess this 
invariance from the very beginning. So we have saved Einstein's equivalence 
principle for our modified model. 

2.3 Some Comments 

• The Hamiltonian H is according to ( |32|) given by 

H = / £x ( 9 ^( dkU f - 9i) + Odifaun)) . (36) 

On the subspace of static solutions of the EOM we have from (26), 
resp. (34) 

diO = 0, resp. U^uf = 9i (37) 

and so we obtain 

H = 8^G J d " Xg ^ > °' (38) 
i.e. the energy of the darkon fluid is positive definitive. But for the 
time-dependent solutions we can always find initial conditions that cor- 
respond to H being negative. However, as we argued in (P, section 
IV) this does not lead to any instability. 

• The EOM (1331) is a continuity equation with a source term for 
d&i = -A<j>. 

Its important property, when compared with the standard nonrelativis- 
tic Poisson equation for (j), is its dynamical nature ie the time develop- 
ment of A(j) is determined by its initial data. 

• Note that the EOM do not contain Newton's gravitational constant. 
The Lagrangian contains it only as a common factor. Moreover, our 
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modified Lagrangian (13"2"j) does not split into two parts involving a free 
and an interacting term. Only the EOM (1201) (resp. (134"}) ) show such 
a splitting. Clearly this strange property is a consequence of phase- 
space reduction due to ( 1T91 . We know only one other physical system 
possessing this property: the interacting Chaplygin gas if one eliminates 
the particle density from the Lagrangian (see [5], section 2.1, item(i)). 



3 Zero mass Galilean symmetry 

The "darkons" introduced in [1] are free massless Galilean particles. As such 
they are a dynamical realisation of the unextended D = (3+1) Galilei algebra 
(for details see [2]). This time the situation is different. The easiest way to 
observe this is to go to the Lagrange formulation. Then using the equations 
derived before we note that the Lagrangian is given by 

n o f (9 



AttG 



d 3 !i (|) m,t),t) -&(a?(|;t),t)) - (39) 



d 3 xg^(x, t). 
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This shows that L does not possess a particle limit. Darkons exist only 
as a self-gravitating fluid. The questions therefore arises which symmetry 
algebra is dynamically realised by such a fluid. 

In this section we are going to consider the Galilei algebra leaving the 
possible conformal extension to the next section. So, let us start with the 
Hamiltonian H (generator of time translations) which has already been given 

by m- 

We can now construct the conserved generators for space translations and 
Galilean boosts. According to Noether's theorem we obtain for the generator 
Pi of space translations 

Pi = d 3 xnpi, (40) 



where 

(npi)(x,t) = -^^(d k g k diU + ndi9) (41) 
11 



and the Galilean boosts are generated by 

K t = tPi + ^ J d 3 x Xl d k g k . (42) 

By using the EOMs it is easily seen that both Pj and Ki are conserved. 
However, the Poisson bracket of these generators vanishes 

{Ki, P 3 } = (43) 

as is easily inferred from the property of Pj as the generator of space trans- 
lations 

{P, A(x,t)} = d t A(x,t) (44) 

for any fluid field A. 

The vanishing of the right hand side of fj4*3|) tells us that the total mass 
of our darkon fluid is zero. 

If we consider, in addition, the conserved angular momentum (generator 
of space rotations) 

J = J d 3 xx x np (45) 

we clearly see that the ten generators H, P i} Ki and Jj satisfy the unextended, 
i.e. zero-mass Galilei algebra. 

4 Anisotropic scaling 

Let us note that we have as a conserved dilation generator 

1 1 f 3 f 

D = Ht - / d s xudigi - - / d 3 xxiupi, (46) 

which leads to the Poisson bracket with the translation generator satisfying 
{Pi,D} = |Pj and thus corresponds to the dynamical exponent z = | ie 
time and space coordinates scale as (see [2]) 

x* = Xxi, t* = XH. (47) 
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To work out the scale dimension z& of a field A(x, t) E (n, u, (f>, 9) defined by 

A*(x,t) = \ ZA A{x*,t*), (48) 

we look at the scale invariance of the action. 

From the requirement that all terms in the Lagrangian scale in the same 
way we find 

z u = z — 2, z^ = 2z — 2, and zg + z n = 3z — 2. (49) 
However, the measure of the integration in the action scales as X~ 3 ~ z and so 

5 
3 



again z = \ and therefore we find that 



1 4 

z u = ——, z <p = zq + z n = 3. (50) 

Note that fl50|) leaves z n undetermined because D is determined only modulo 
a term ~ Q, where Q is the conserved charge Q = J d 3 x n6. In (|46p we have 
fixed this term in such a way that we have scale invariance of the particle 
number 

N = J d 3 xn(x,t) (51) 

and so we get z n = 3 and in consequence zg = 0. We conclude that the 
symmetry algebra of our darkon fluid becomes the expansion-less conformal 
Galilean-type algebra with z = | (cf. [2]). 

To derive the form of our EOM , either for the scale invariant solutions or 
for cosmological perturbations, it is convenient to perform the transformation 



(x,t) ->■ (w,t) (52) 



with w = x ( a and let 

_ 3z^4 

A(x,t) = (J-J ' A(w,t). (53) 

As w is scale invariant we see from ( 148)) that A scales as 

A*(w,t) = A{w,\h) (54) 
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and so A is scale invariant (A = A*) if 



d t A(w,t) = 0. (55) 



But, due to (1531) . the partial derivatives of A transform as 



t N -!(i+* A ) 



and 



$A(£,t) = ( - ) ^A^t) (56) 



<M(^)=(£) 5 U- ^K9 A + ^) ) A(w,t). (57) 



Hence, in terms of A the EOM (13"4T), (Ell). (j2Bjl and (1551) take the form 

3 1 t 

(t<9 4 - -w fc d fe + -)tt + -^-{d k u) 2 = -to4> (58) 



and 



td t n + d k ({t d k u - ^w k )hj = 0, (59) 

td t 9 + (t d k u - ~w k ) d k 6 = (60) 
5 



(td t - \w k d k - 2)A0 + t di{&4>diu) = t o di(hdi0). (61) 
5 



It is striking that these EOM contain an explicit time dependence in 
the form of the term td t . To give them the standard 1-st order form it is 
convenient to perform the transformation 

r = log (f) , (62) 



to 

or, for the derivatives, td t = d T . 

Then the EOM ( 1581I6T!) are derivable from the Lagrangian (note that 
and not (f> is the relevant variable in L) 

If/ 3 2 

L = £ / d 3w [di(d T + (t d k u - -w k )d k - - (63) 

~-^9i + 0(d T n + d k ((t d k u - ^w k )n] 
14 



which may be obtained from the original action by inserting the transfor- 
mation ( )52ll53l) and ( 1621) . The corresponding Hamiltonian H can be read off 
from ( I63p by writing it as 

L = l — [ d 3 w(gid T diU + 6d T n) - H. (64) 

47rG J 

Now we have to be cautious: The Hamiltonian H is not the transformed form 
of H (eq. ( 156"|) ) but it comes from inserting the above mentioned transforma- 
tions into the dilation generator D (|46|) . The reason for this replacement of 
H by D is an immediate consequence of the transformation (IBIZl) which leads 
in Q5U) to 

A*{w,t) = A(w,r+ - log X) (65) 

3 

i.e. time rescalings become time translations (both to be taken at fixed w). 

Note that for scale invariant solutions, characterised by (1551) eq. (1601) 
turns into 

tod k u - -w k ) d k 9 = 0. (66) 
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Thus we have to distinguish three different possibilities 
• 1. 

3 

«o(0 = Y^o 1 ^ 2 ' ( 67 ) 

where w = \w\. Then f)59p is identically fulfilled and from (158 p . resp. 
(IBTj) we obtain 



^tfw\ (68) 



and, resp., 

18 

d k (n (w)d k 6 (w)) = "Y^^o 3 - (69) 

In the next section we shall show that the choice ho(w) =const. leads 
to the cosmological solutions for the early universe. 



d k 9 = (70) 
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but t Vu — ~w ^ 0. Then ( 1591) has no radially symmetric solutions 
valid for all w > 0. To prove this statement we note that the equation 
dk(wkf(w)) = has, for nonvanishing w, the unique solution }'{w) ~ 
w~ 3 leading to dk(wkf(w)) ~ 5(w). We have not so far succeeded in 
finding, analytically, the anisotropic solutions. 

3- 

t Vu - ^w^j JL V0. (71) 
The discussion of ( 17T|) lies outside the scope of the present paper. 

5 Cosmological solutions and their perturba- 
tions at early times 

The cosmologically relevant solutions, as derived and discussed in p], are 
given by 

M*,*) = ^fy (72) 
uo(x,t) = -x, (73) 

and 

3DG 

g {x,t) = —g(a)x, (74) 

or 

where D is a positive constant. 

The unknown cosmic scale factor a(t) and the function g(a) obey the 
following set of differential equations: 

a = -^-g(a) (75) 



a 2 



and 



9(a(t)) = 4 (76) 
or 



where /? is a constant. Then we have 



y.2 



a 2 



0„(x,t) = 2vrG/3-. (77) 
16 



The last expression (177)) comes from the observation that for the cosmo- 
logical solutions the equation for 9 reduces to 

d t + -x k d k )e = 0, (78) 
a J 

which has (1771) as its solution. The overall constant comes from inserting the 
expressions above into the equation for g^ ie (1351) . 

Comparing (I72H77I) with the scale invariant solutions of case 1. ( 16711691 
we note that we have agreement if 

3-D 

no(w) = — , (79) 



and 



a(t) = ( 1 V , (80) 



where to becomes the cosmological time at present and 

According to pQ eq. (180)) gives the cosmic scale factor a(t) at early times, 
i. e. our cosmological solutions are scale invariant for times before structure 
formation sets in. Therefore the scaling law (14 7p is not applicable to the 
relations between microscopic scales and galactic ones (see [TO]). 



5.1 Perturbations 

Given the cosmological solutions we can now consider their perturbations. 
To do this we take any fluid field A 6 (n, U{, gi, 9) and put 

A(w,t) = A (w) + 8A(w,t) (82) 

where the first term in (182)) denotes the cosmological background solution and 
the second its perturbation. Note that 8n is proportional to the fractional 
density perturbation 5h ~ ^ and w is the coordinate comoving with the 
background. 
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Putting the expressions (182]) into the Lagrangian (I6"3l) and keeping terms 
quadratic in perturbations gives us (we abreviate A = 5 A) 

aIg j ^ { m + ~ 2 # + I^ 1( ^ )2 + (83) 



A 3D A . 87r „ 
+ 9{d T n + — t A«) + i25Dtl nW U 

It turns out that this Lagrangian is invariant with respect to the rescal- 
ing of the comoving space coordinates w when the fields A G (n,9,gi,u) 
transform as 

A(w,t) A*{w,t) = A r?A i(Aw,r) (84) 

with 

1 3 5,1 

Vu = 2' ^ = 2' ^ n = 2 7/9 = 2" ^ 

By Noether's theorem we obtain for the generator .D s of this symmetry 
D s = J d 3 w \ -digi(~ + w k d k )u + + w k d k )nj . (86) 
The EOM following from the Lagrangian (183"]) are 

(Or + g)« = -*o0, (87) 
3D 

d r n + — t Aw = 0, (88) 

47T 



Sn 

9r0 = U5Dt 2 W k d kU, (89) 



and 



( * - l )A * + 2^ Afi - ^<» A * + i2^ a - (tt, - fl) = °- (90) 

Applying d T to ( 190]) it is possible to eliminate n and 6 1 . One then obtains 

1 18 12 

<9 T (d T - -)A0 - — A0 - -—(4 + ^)A« = 0. (91) 

With the ansatz 

u(w,t) = e aT u x (w), (92) 
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for some a we find from (19 ip with (I87p that 



d 

(ai + a 2 w—)Au 1 (w) = 0, (93) 

where the constants are given by 

/ 2 19x 30 12 
ai = a{a ~ 25 } + 125' a2 = 125 (94) 

Let us now consider the particular case when u\ is invariant with respect 
to the space rescalings Then we have in spherical coordinates 

ux{w) = w-$f(9,<p), (95) 

which gives us from (|93l) an equation for a 

a x - ^a 2 = (96) 



and therefore 



«(« 2 = (97) 



with the solutions 

fig 

a = 0, a = ±W— (98) 
V 

Expanding f(8, (p) into spherical harmonics Y^ m we obtain for the growing 
mode of h 

/ , \ 0.8718 

n(^,t) ~ f-J w;-3^ >m . (99) 

As the EOM (I87H90P are linear the general solution for h is a linear com- 
bination of the partial solutions (19"9"|) with their relative weight determined 
by the initial conditions. 

Clearly, the existence of structures, such as galaxies and stars, in the 
universe suggests that the perturbations, at early times, were not radially 
symmetric (cf |llj). 
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6 Adding higher-order derivative terms 

A generalization of the EOM (|20|) which preserves gauge invariance is ob- 
tained by adding terms containing spatial derivatives of gi. Such terms are 
generated by adding to the Lagrangian (132|) terms of the form (we consider 
only terms which are quadratic in g^) 



The first trem in L' corresponds, for = —8$, to the space-part of the 
interaction term for a Lifshitz scalar field (see [6] and the literature cited 
therein). Clearly [K n ] = (length) 2 ". 

Note that the EOM remain the same as before except for the Euler equa- 
tion which now becomes 



Note that 

• The terms proportional to the K n do not contribute to the cosmological 
solution given by g k (x,t) = Xkg(t). 

• Whether the K n vanish or not the system of equations does not possess 
static radially symmetric solutions. 

• Note that fllOOH scales differently than the gf term in (|32|) . To be 
specific suppose that we consider only the first term in fllOOp as the 
interaction term. Then we find that our fields have the scale dimensions 
z u = |, z g = | and the dynamical exponent is given by z = |. This 
suggests that when we have both interaction terms we have a model 
which has a dynamical exponent z = ~ at short distances and z = | at 
long ones. 

7 Conclusion and Outlook 

The main aims of this paper involved: 




(100) 



di(d t u + -{d k uf) = g { - K^d^k) + K 2 diAd k g k + .... (101) 
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• the elimination of the unphysical degrees of freedom in the original 
formulation of the model pQ for the description of a self-gravitating 
darkon fluid, 

• the demonstration that this darkon fluid exhibits the zero mass Galilean 
symmetry and, in consequence, possesses anisotropic scaling symmetry, 

• the exploitation of the anisotropic scaling symmetry for the formulation 
of cosmological perturbations within the dark sector of the universe. 

In order to obtain results that could be compared with cosmological data 
we have to add baryonic matter and the higher spatial-derivative terms as 
further perturbations. Then, however, everything becomes more complicated 
as any of these contributions break the dilatational symmetry of our EOM. 
Thus we will leave such an extension of our present work to a subsequent 
paper. 

Our model possesses the symmetry with respect to time-dependent space 
translations 

Xi -> x'i = Xi + ai(t). (102) 

We plan to try to generalize our model so that it comes as close as possible 
to General Relativity which possesses symmetry with respect to arbitrary 
space-time diffeomorphisms (Diff). In a nonrelativistic setting we have 
therefore to consider the subgroup of Diff which leaves the Newtonian time 
fixed (modulo reparametrizations). 

Xi x[ = Xi + a,i(x,t), t — > t' = f(t). (103) 

In order to do this we have to start by 'covariantizing' the time derivative 
terms in our Lagrangian fl32|) . Such a 'covariantizing' procedure has been 
worked out already some time ago for massive point particles in (1+1) [12J, 
respectively (2+1) (13] dimensions. This project is currently being studied 
and we hope to present our results soon. 
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